Minigap in a SN junction with paramagnetic impurities 
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We study the effect of spin-flip scattering on the density of states in a long diffusive S-N junction 
with transparent interface. We calculate the critical value of the spin-flip scattering rate at which 
the minigap closes and give the dependence of the minigap on the spin-flip rate. For the system we 
consider, the minigap and the critical spin-flip rate have the scale of the Thouless energy, and not 
of the superconducting gap. 



Proximity structures made of superconducting and 
non-superconducting elements in contact with each other 
is a topic of current theoretical and experimental inter- 
est. Most of the experimental applications concern dif- 
fusive materials for which an appropriate formalism has 
been derived^ In the diffusive regime the motion of the 
electrons is governed by frequent scattering on impurity 
atoms: the elastic mean free path l e is much smaller 
than the relevant length scales of the system, namely the 
superconducting coherence length £ and the dimension 
of the normal (N) and superconducting (S) electrodes. 
As a consequence of Andreev reflection and its interplay 
with disorder, a gap is opened in the spectrum of the 
N-region (see for example 2 *^*^ and references therein). 
This minigap does not depend on the position in the 
junction and is of the order of magnitude of the Thouless 
energy i?Th = ^s, where D — ^^y^ is the diffusion con- 
stant and L is the length of the normal wire. Here and 
in the following, we use the units with % = 1. 

Elastic scattering on impurities is not harmful for 
isotropic superconductivity and is directly connected to 
the appearance of the minigap. On the contrary, scat- 
tering on magnetic impurities leads to the destruction 
of singlet superconductivity. The exchange interaction 
between the spin of the electrons and the magnetic im- 
purities breaks the time reversal symmetry between the 
electrons in Cooper pairs and reduces the superconduct- 
ing correlations. Scattering on magnetic impurities can 
be described in terms of a spin- flip scattering rate T s f. 
In the Born approximation^! T s f is proportional to the 
impurity concentration. The energy gap E g in the elec- 
tron density of states (DoS) of a bulk superconductor is 
lowered by spin-flip scattering 



and closes for the critical concentration of magnetic 
impurities^ Fg f ulk = -y- Here, the order parameter A 
itself depends on T s f. 

In the "dirty" limit, the equations of motion for the 
quasiclassical matrix Green function (averaged over the 
fast Fermi oscillations and the momentum directions) can 



be reduced to the Usadel equations^ Usadel equations 
have been very useful in the analysis of position depen- 
dent problems like proximity effect. Spin- flip and inelas- 
tic scattering processes can be included in the Usadel 
equations through the corresponding scattering rates T s f , 
respectively Tiniifi 

The impact of spin-flip scattering on the proximity 
effect and the minigap can be controlled by selectively 
doping the normal and/or superconducting regions with 
magnetic impurities or applying a magnetic field to a thin 
film or wire^ii Belzig et aii^ have numerically shown 
that the minigap is reduced and vanishes for large values 
of r s f . For the case of a junction with L = 1.1 £ (corre- 
sponding to the Thouless energy of the order of A), they 
find that the minigap closes when T s f as OTA. Other 
worksi&i^ have analyzed the effect of magnetic impuri- 
ties in the transport properties of S-N junctions, finding 
that the Thouless energy is the scale which controls the 
effect of spin flip scattering on the proximity effect. Due 
to the self-consistent treatment and the non-linearity of 
the Usadel equations the complete numerical solution of 
the problem is difficult. It is important to obtain sim- 
ple expressions which allow comparison with experiment 
without solving the full equations for each set of param- 
eters. 

In this note we obtain such expressions for the case 
of a long junction. We show that a finite, but small T s f 
suppresses linearly the minigap which closes at « 
0.62 i?Th- We also present an asymptotic expression for 
the dependence of the minigap on the spin-flip rate near 
the critical value. Above this value, we calculate the zero 
energy density of states. 

We consider a finite size normal metal N, of length L 
connected to a semi-infinite superconducting terminal S 
by a transparent interface. We assume that electronic 
motion is diffusive in both the normal and supercon- 
ducting parts. We restrict our discussion to a quasi one- 
dimensional geometry /wire or to a thin film and neglect 
any dependence on transversal coordinates. The origin of 
the coordinate x is fixed at the S-N interface. Introduc- 
ing the angular parametrization 2 g R = cos 9, f = sin# 
for the normal, respectively the anomalous component of 
the retarded Green function, the Usadel equations in the 
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N-region, where the pairing interaction vanishes, become 
Ddl a 



{ie - r in ) sin 9 - 2T sf cos 9 sin 9 = 0. (2) 

The proximity angle 9 is a function of the energy e and 
the position x. The electron DoS, in units of the normal 
state bulk value v$ = is given by 



Re [cos#(e, x)] 



(3) 



The boundary conditions for the quasiclassical equations 
have been discussed in Refs. Ilallrl At the interface with 
vacuum, the conservation of the quasiparticle current 
yields 



d x 9{x = L) = 0. 



(4) 



We will study these equations analytically using sim- 
plified boundary conditions at the S-N interface, where 
we impose the superconducting bulk value (at zero en- 
ergy) of the proximity angle 



9{x = 0) = 



(5) 



This boundary condition is justified for energies much 
smaller than the superconducting order parameter A and 
if the normal part is much more disordered than the su- 
perconducting part. 5 

Since the scale of the superconducting order parameter 
A does not appear in the rigid boundary condition JSJ), 
we will write in the following the energies and the length 
in units of the only other relevant scale for our system: 
the Thouless energy £/rh, respectively the width of the 
N-region L. 

The boundary conditions (1 1I5[I and the calculations 
presented here for the S-N junction can also be applied 
to describe S-N-S junctions with no phase difference 
between the superconducting terminals^ It is important 
to pay attention to the choice of the energy scale i?xh- 
A S-N-S junction of length unity is equivalent to a S-N 
junction of length i. Therefore all the energies must be 
multiplied by a factor four if we consider S-N-S junctions. 

In the limit of small energies and scattering rates the 
Usadel equation @ can be linearized in terms of the 
deviation 69 from the value of the proximity angle at the 
S-N interface (JSJ. At this level of approximation the DoS 
in the N-region is given by 



T<gx{2-x), 



(6) 



where rjf = r in (l + 2r sf ) and x € [0, 1]. We see that if 
r in = a small spin-flip scattering does not create states 
at small energies. If, on the contrary, Ti n is finite, T s f 
modifies the density of states (constant at small energies) 
induced by the inelastic scattering. 




FIG. 1: Numerical gap curve (solid line): comparison with 
asymptotic expressions (dashed lines). 



To study the presence of solutions where the electronic 
DoS © vanishes we focus in the rest of the paper on the 
situation with Tin = 0. We introduce a new notation 
for the proximity angle 9 below the minigap and write 
9 — j +i(3 with p real. The minigap E g is, by definition, 
the maximal energy compatible with a real (3 and can be 
obtainedi£ii£ using a first integral of J3J 



d x p = 2V/(/3 1 )-/(/3), 



(7) 



where the superscript 1 denotes the value at x = 1 and 
f(P) = t sinh (3 + r s f sinh 2 /?. Integrating equation Q 
over the junction, we get 



1 



dp 



o 2^/(/3 1 )-/(/3) 



f. 



(8) 



Without spin-flip, we recover with this relation the well- 
known valued of the minigap E® ss 0.78. The critical 
value of the spin-flip at which the minigap in the DoS 
closes is (see the discussion in Appendix) 



isf " 16 



0.62. 



(9) 



For energies e < and r s f = equation JSJ) is solved 
for two different values of /3 1 . One of these solutions leads 
to a diverging /3 1 when the energy goes to zero and we re- 
ject it using a continuity argument. This continuity argu- 
ment is commonly accepted in the quasiclassical approxi- 
mation, but the diverging branch may play an important 
role in the discussion of the presence of a non-zero subgap 
DoS resulting from mesoscopic fluctuations^ Consider- 
ing a finite spin- flip scattering rate T s f , we find that the 
second branch of the solution no longer diverges at zero 
energy. If we increase the spin-flip rate up to a critical 
value the two zero energy solutions merge at T^ f . The 
critical value can therefore be determined taking the limit 
f3 1 — > of the integral in the l.h.s. of equation JSJ) at zero 
energy. 

The complete dependence of the minigap on the spin- 
flip rate (Fig. ^) can be obtained by a simple numerical 
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integration of equation JSJ). But it is possible to derive 
the asymptotic form of the gap curve for T a i — * and for 

r S f -> r° f . 

In the limit of a small spin-flip rate, we can expand the 
integrand in JSJl in the small parameter Denoting 
j3 1 ~ 1.421 the value of the proximity angle correspond- 
ing to the zero spin- flip value of the minigap Eg, we find 
the resulting correction 



{E° g - dr sf ) 



(10) 



where the coefficient C\ is given by 



d = 



r/3 1 (sinh^M-sinhgj 
^° (sinh/Si-sinh/S) 172 

Jo 



rig 



3.09. 



(11) 



(sinh/fi-sinh/S) 



The minigap decreases linearly with increasing spin-flip 
scattering rate. From the magnitude of C%, we can see 
that even a small spin-flip rate strongly affects E g . 

To obtain the analytic behavior close to r^ f is more 
tricky as all <5r s f = — T s f, (3 and ^ are small. Fol- 
lowing the procedure detailed in the Appendix we obtain 
the asymptotic dependence of the minigap near the clos- 
ing point 



E„ 



2 5T si 



3/2 



(12) 



In Fig. ^ we compare the asymptotics l|10[) and H12|) 
with the numerical gap curve. Dimensions are reintro- 
duced in the graphs for clarity. 

For r s f > r^f the DoS in the N-region is finite at any 
energy. In this domain the Usadel equation at zero 
energy has a solution with 9 real. Applying again the 
procedure of the first integral (JJJ , but this time for a real 
9, we get 



8(x) 



(19 



Vcos 2 9 1 — cos 2 9 



= -2Vr sf x 



(13) 



Inverting this elliptic integral, we can find a complete zero 
energy solution for equation J2J). This solution involves 
the function dn(u, k), which is one of the Jacobi elliptic 
functions, defined as inversions of the canonical forms of 
elliptic integrals (we follow the notations used in Rcfs. 5, 
2C 



suit 



dn(2y/T7f(x- l),cos9 1 ) 



(14) 



where 9 1 can be determined imposing the rigid boundary 
condition (J5J at x = 



sinfl^dn -2v/frf , cos 9 1 



(15) 
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FIG. 2: Dependence of the zero energy local DoS, at the open 
boundary, on the spin-flip scattering rate: asymptotic expres- 
sion near (solid line), complete elliptic solution (dashed 
line) saturating to the normal state bulk value vq and its 
asymptotics at large T a f (dots). The inset shows the depen- 
dence on position of the DoS for F s f = Et\i and T s f = 5i?Th- 



In the inset of Fig. [21 we use relation (|14|) to represent 
the dependence on position of the DoS in the N-region for 
two different spin- flip rates. As expected the density of 
states increases when we move away from the interface. 
For large spin- flip scattering rates, the elliptic solution 
(|14fl saturates to the normal state bulk value 9{x) = 0, 
everywhere except in a thin domain close to the S-N in- 
terface (rigid boundary at x = 0). 

Near the critical spin-flip rate, we found a square root 
dependence of the zero energy local DoS on r s f . Expand- 
ing the integrand in (|13|l in the small cos 9 1 , we obtain 
that for r s f — ► r^f the density of states at the interface 
with vacuum is given by 



v(e=0,x=l) 



V2 



5T s t 



1/2 



11 

8V2 



ST, 



s[ 



■ sf 



3/2 



. (16) 

In the limit of large spin- flip rates, when the density of 
states approaches the normal state one, the expansion 21 
of the elliptic integral in l|13fl near cos 9 1 = 1 leads to an 
asymptotic expression for the DoS at the interface with 
vacuum 



!^(e=0,a;=l) 



1 



-4v/rv 



(17) 



In Fig. |3 we compare the expressions (|16I17|I for the 

DoS at the open boundary with values obtained using 
the complete zero energy elliptic solution (|14l) . 

In summary, we have studied the effect of spin-flip scat- 
tering on the minigap and zero energy density of states 
of a normal slab connected to a superconducting elec- 
trode. We have obtained analytic expressions for the 
dependence of the minigap on the scattering rate both 
for small spin-flip scattering rates and close to the crit- 
ical value at which the gap closes. This critical value is 
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controlled by the Thouless energy. For values larger than 
r° f we find finite DoS at zero energy with a square root 
dependence on T s f. These results are valid for long S-N 
or S-N-S junctions with transparent interfaces, at small 
energies. The results on the minigap are not restricted to 
the case in which there is no spin-flip scattering in the su- 
perconducting regions. The only requirement is that the 
Thouless energy of the normal part is small compared to 
the gap in the superconducting one. These results allow 
easy comparison with experimental measurements. 

We thank M. Sigrist for useful discussions and acknowl- 
edge financial support from NCCR MaNEP of the Swiss 
Nationalfonds and the Spanish Science and Education 
Ministry through Ramon y Cajal contract. 



find the maximum value of Y over z\ for a given value of 
the parameter a. 

The critical spin-flip scattering rate © can be ob- 
tained setting a — and maximizing (|A.1J| . It turns 
out that Y(zx, a)| Q =o is largest for z\ = 0. The next step 
is to go to finite energies and expand Y in a. We write 



dY{z\ 1 a) 

Y{zi,a) « Y(zi,a)\ a=0 H 5 \ a =o<x> 



da 



where 



dY( Zl ,a) 
da 



1 



Zl 



l + z{ 



(A.2) 



(A.3) 



APPENDIX: GAP CURVE CLOSE TO T= f 

In this appendix, we derive the asymptotic form of the 
minigap curve close to r^ f . We have seen previously that 
the minigap E g is the largest energy compatible with 
equation JHJ. Introducing 



z 
z\ 

with a 



sinh ft + a 
sinh (3 1 + a 



2IV 



we can rewrite equation JHJ in the form 
1 1 



2vr sf = 



y/l + (z - a)- 



The integral in the r.h.s. is a function of z\ and a, which 
we denote Y(z\,a). To determine the minigap, we will 



For a small a, the maximum of Y is expected to be close 
to the zero energy value Z\ = and the second term 
~ O(zi) can be neglected. The first term in the r.h.s. 
can be expanded 



Y(zi,a)\ a =o 



1 



o VT 



(z lS f 



ds. (A.4) 



where s = . Taking the derivative of Y over Z\ , using 
the expansion l|A.2|l , we find that the maximum of Y is 
obtained for 



1/3 



Zl 



(A.5) 



-dz (A 1) n I 

\ ■ J substituting back this result in (|A. 2|) and using the defini- 



tions of a 



2F 3{ 



and Y 



/ T s f we obtain an expression 



for E g (T s {). Finally we write T s f = T% — 5T S { and expand 
in the small 5T S { to get the asymptotic dependence <|12ll • 
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